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Abstract 

A numerical study was performed to simulate the critical 
flow through a supersonic inlet. This flow field has many 
phenomena such as shock waves, strong viscous effects, tur- 
bulent boundary layer development, boundary layer separa- 
tions and mass flow suction through the walls, (bleed). The 
computational tools used in this study were two full Navier- 
Stokes (FNS) codes.. The supersonic inlet that was analyzed 
in this study is the Variable Diameter Centerbody, (VDC), 
inlet. This inlet is a candidate concept for the next genera- 
tion supersonic transport. Application of the code to the 
inlet geometry involved effort in generating an efficient grid 
geometry and specifying boundary conditions, particularly 
in the bleed region and at the outflow boundary.. Results for 
a critical inlet operation compare favorably to Method of 
Characteristics predictions and experimental data. 


Nomenclature 

Principal Symbols: 

Unless specified otherwise, all variables are non-dimension- 
al ized. 

A flow area, m 2 (ft. 2 ) 

A c capture area, 

0.7073 m 2 (7.61 ft. 2 ) TF-30 size 
0.1400 m 2 (1.51 ft. 2 ) J-85 size 
c p specific heat at constant pressure 
d height from surface to center of 
total pressure tube, m (ft.) 

E t total energy per unit volume, 

k effective thermal conductivity coefficient 

H annulus height at local diffuser station, m (ft.) 

h distance from surface, m (ft.) 

L length, m(ft.) 

M Mach number 

m/m c mass flow ratio 


p 

pressure, N/m 2 (psf) 

p 

static pressure 

Pr 

Prandtl number 

R 

Gas constant 

R c 

inlet capture radius, 

0.4745 m (18.68 in.) TF-30 size 
0.21 16 m (8.331 in.) J-85 size 

r 

radius, m (ft.) 

Re 

Reynolds number 

t 

time 

T 

Temperature, K(°R) 

u,v,w 

Velocities in the x, r and swirl directions 

X 

axial direction, m (ft.) 

X 

axial direction 

Y 

ratio of specific heats 

0 

circumferential angle 

X 

second coefficient of viscosity 

V 

viscosity coefficient 

P 

static density 

T 

shear stress 

Subscripts: 

bl 

centerbody bleed 

by 

bypass 

c 

capture 

cf, ex 

diffuser exit or 
compressor face station 

inf, oo 

freestream 

lam 

laminar 

max 

maximum 

min 

minimum 

ref 

reference or normalizing conditions 

turb 

turbulent 

X 

local 
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This study’s goal is to accurately predict critical inlet opera- 
tion with an advanced computational tool. The main chal- 
lenge to accomplish this goal is the accurate modeling of the 
inlet’s throat bleed region. Additionally this work demon- 


strates the capabilities of a new ‘general purpose’ fluid flow 
solver, PROTEUS, for use in analyzing inlet performance. 

The inlet of a supersonic cruise aircraft is a crucial element 
in an efficient propulsion system. The inlet compresses the 
flow entering the engine with shock waves. For good perfor- 
mance, a supersonic cruise inlet is designed to compress the 
freestream flow to a nearly sonic Mach number, (Ml. 2 to 
Ml. 3), at the inlet’s minimum cross-sectional area or throat. 
Inherent in this compression process are disturbances in the 
airflow along the walls of the inlet, a problem known as 
shock wave/boundary layer interaction. Due to the critical 
Mach number in the throat, even a slight disturbance can 
radically alter the inlet flow, i.e. inlet ‘unstart’ or even 
‘buzz’. ‘Unstart’ or ‘buzz’ are inlet conditions in which the 
normal shock in the throat has been expelled forward in 
front of the inlet and becomes highly unsteady. To keep the 
flow well conditioned at critical inlet operation, the dis- 
turbed airflow is removed through small openings in the 
inlet walls in a process known as ‘bleed’. In short, an effi- 
cient supersonic inlet could not operate without bleed. Due 
to the complexity of the airflow, design of inlet ‘bleed’ is 
largely based on empirical guidelines. 

To develop a better understanding of the flow complexities, 
a detailed test program and a computational solution to the 
complete fluid flow equations is needed. In fact, this study is 
specifically targeted at the Variable Diameter Centerbody, 
(VDC), Inlet which is a design for which experimental data 
exists and future wind tunnel testing is planned. The VDC 
inlet was designed for a M2.5 cruise passenger transport. 
Computational analysis tools used to prepare for wind tun- 
nel testing has and will impact the understanding the inlet 
flowfield. 

The application of Computational Fluid Dynamics, (CFD), 
to the problems in supersonic inlet flow was begun with 
Euler solvers over twenty years ago. The inlet in the present 
study was in fact designed with the aid of a Method of Char- 
acteristics code, reference 5. Further progress in applying 
CFD to the strongly viscous flow field in the throat and dif- 
fuser of supersonic inlets was delayed until the development 
of full Navier-Stokes, (FNS) codes and much faster comput- 
ers became available. Some successful studies have been 
made in this area but are limited because they require exper- 
imentally derived knowledge of the bleed flow rates and/or 
pressure distribution in the bleed area. These efforts take a 
traditional approach to modeling the bleed flow by setting 
boundary conditions at the wall that allow flow to exit the 
computational domain, references 12 to 17. 

Two studies have taken a different approach by actually 
solving the flow from the throat through the wall boundary 
layer into the bleed slot, references 18 and 19. This 
approach is very time-consuming but has shown some 
insight into the bleed flow region. The conclusion is that 
bleed cannot be modeled by simple boundary conditions. 
The insight from this work can be used to devise more accu- 
rate albeit complicated boundary conditions. This insight 
will be used to guide bleed modeling in this study. 


Beyond the scope of this work, the ultimate goal is to use 
CFD as a tool to investigate the ‘goodness’ of various inlet 
designs. With several orders of magnitude increase in com- 
puter capability or algorithm speed, (in perhaps 5 years - 
maybe 507), inlet design may be accomplished largely on a 
computer with wind tunnel testing used only for verification 
purposes. 

II. Background 

A. Variable Diameter Centerbody 

The Variable Diameter Centerbody, (VDC), inlet is an 
ongoing research program at the Lewis Research Center. 
The VDC inlet is a mixed compression, axisymmetric inlet 
that has potential application on the next generation super- 
sonic transport. This inlet was identified as one of the most 
promising axisymmetric concepts for supersonic cruise air- 
craft during the SCAR program in the late 1970’s, reference 
7. Some of its features include high recovery, low bleed, 
good angle-of-attack tolerance and excellent engine airflow 
matching. These features have been demonstrated at Lewis 
in the past by the design and testing of fixed hardware mod- 
els, references 1 to 4. A current test program in the LeRC 
10’xlO’ Supersonic Wind Tunnel, (SWT), will attempt to 
duplicate these features on flight-like hardware that actually 
incorporates the variable geometry. 

An outcome of the supersonic cruise research (SCR) pro- 
gram identified the VDC inlet as an important technology 
thrust to continue funding, reference 6. It is an axisymmet- 
ric inlet with 45% internal contraction ratio needed for high 
performance at its design cruise Mach number of 2.5, figure 
1. Other aspects of the inlet design include the variable 



diameter centerbody and a focussed cowl compression on a 
slotted bleed region in the centerbody. The variable diame- 
ter centerbody allows large variations in throat area and air- 
flow to provide good compatibility with the engine. The 
focussed cowl compression minimizes bleed flow require- 
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ments and reduces the inlet length and resulting weight. 
Aerodynamic testing of the concept was done with fixed 
hardware in the early 1970’s and verified the high expected 
performance of this concept. An schematic view of this 
model is shown in figure 2. The model had centerbody and 



FIGURE 2., Fixed-hardware model, 
schematic view 

cowl bleed for performance and shock stability and over- 
board bypass air for engine matching. Vortex generators 
were installed downstream of the throat to prevent separa- 
tion in the subsonic diffuser. This model was tested in the 
LeRC 10’xlO’ Supersonic Wind Tunnel at Mach numbers 
of 2.5 and 2.0. For economic reasons the mechanical design 
of that test inlet was simplified to incorporate fixed center- 
body configurations. The fixed-hardware inlet model was 
sized and tested with a TF-30 turbofan engine. The cowl lip 
radius, R c , was 18.68 inches. A photo of the model installed 
in the 10’xlO’ SWT is shown in figure 3 



FIGURE 3., Fixed-hardware model. 

Installed In 10’xlO’ SwT 

The essential features of the inlet design incorporate a 
bicone centerbody of 12.5° and 18.5° half angle cones and 


an initial internal cowl angle of 2°. The design philosophy 
for this mixed compression inlet is to utilize a bicone spike 
to provide the maximum external compression compatible 
with high total pressure recovery and low cowl drag. As a 
result, 45 percent of the supersonic area contraction is inter- 
nal for the Mach 2.5 design condition. 

Mechanical design of the VDC inlet with the variable 
geometry began in 1982 and a complete set of drawings was 
finished in mid- 1984. Unfortunately, programmatic restruc- 
turing canceled the program with only a fraction of the 
hardware fabricated or procured. The High Speed Research 
program has revived interest in a commercial supersonic 
aircraft in general and this inlet program in particular. The 
test program in the LeRC 10’xlO’ SWT is slated to begin in 
the summer of 1992. 

The VDC inlet incorporates an umbrella-like mechanism to 
create a variable diameter centerbody, figure 1. The vari- 
able-hardware model is sized for a J-85 turbojet engine and 
is less than half the size of the fixed-hardware model, 
R c =8.31”. The supersonic diffuser of the J-85 sized VDC 
inlet is geometrically scaled from the fixed hardware model. 
The subsonic diffusers are slightly different. 

This analytical study was undertaken, in part, to prepare for 
the experimental test program. Input for the analysis was 
setup for the fixed-hardware, TF-30 sized, inlet. Compari- 
son between experimental test results from the fixed geome- 
try model and various computational fluid dynamic analyses 
will be made. 

B. PROTEUS Code 

A computer code, called PROTEUS, is used to solve the 
two-dimensional, axisymmetric, Reynolds-averaged, 
steady compressible Navier-Stokes equations for the flow 
through the VDC inlet at its design Mach number of 2.5. 
The PROTEUS code was developed for aerospace propul- 
sion applications and was designed to be easy to use and 
modify. Code readability, modularity, and documentation 
were emphasized in developing PROTEUS, reference 8. 

The governing equations are written in Cartesian coordi- 
nates and transformed into generalized nonorthogonal 
body-fitted coordinates. They are solved for this application 
by marching in time using a fully-coupled altemating- 
direction-implicit solution procedure with generalized first 
order time differencing. The boundary conditions are 
implicit and all terms, including the diffusion terms, are lin- 
earized using second-order Taylor series expansions. Turbu- 
lence is modeled using a Baldwin-Lomax based algebraic 
eddy viscosity model. 

The program has many options which may be selected 
through an input file. For this study both Euler solutions and 
full Navier-Stokes solutions were obtained. The energy 
equation was eliminated through the assumption of constant 
total enthalpy. Explicit nonlinear artificial viscosity was 
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used to damp the pre- and post-shock oscillations in the 
supersonic regions of the flow. The computational grid was 
generated external to the PROTEUS code due to the com- 
plexity of the inlet geometry. 

The PROTEUS code has many options not discussed here 
but fully documented in reference 8. This reference covers 
the analytical development, user application with test cases, 
and detailed program information of the code. One of the 
main capabilities of the code not used in this study is that it 
can be used in a ‘time-accurate’ mode for unsteady flows by 
using second order finite differencing in time. Future capa- 
bilities in progress are a three-dimensional version of the 
code. Also in the planning stages are a multi-blocked ver- 
sion allowing a flow field to be broken up into many grids. 
This allows problems requiring a large number of grid 
points to be solved as a series of smaller interacting prob- 
lems on computers with finite memory resources. Blocking 
also relieves the waste of gridding non-flow areas such as 
struts. 

C. PARC Code 

The PARC code is a well established full Navier-Stokes 
code originally developed at Arnold Engineering Develop- 
ment Center as a general purpose robust code to reduce the 
costs and risks in test planning, execution and analysis, ref- 
erence 9. Subsequent use has included inlet and nozzle anal- 
yses done at NASA-Lewis. In many respects, the code is 
similar to PROTEUS so that much of the initial setup could 
be directly transferred to the PARC code analysis. In fact, 
the PARC code analysis was begun using solution and grid 
files from the PROTEUS analysis. As such, only the major 
differences between the codes will be pointed out here. 

The algorithms within the flow solvers of the two codes dif- 
fer somewhat. The PROTEUS code essentially uses the 
Beam -Warming formulation that must solve a blocked pen- 
tadiagonal set of matrices to advance the solution in time. 
Using Pulliam’s approach, the PARC algorithm begins with 
this setup, factors and uses LU, (Lower-Upper), decomposi- 
tion on the matrices to reduce the complexity of the solu- 
tion. This approach is inherently faster than the algorithm 
used by the PROTEUS code and offers a factor of 2 increase 
in speed for the present application. In addition, the PARC 
coding is tailored for the CRAY supercomputer and takes 
six times less memory to run the present application. While 
both codes have options to provide ‘general’ boundary con- 
ditions, the PARC input is tailored to set large regions of the 
boundary to the same condition while PROTEUS allows 
grid point by point specification of the individual thermody- 
namic quantities. This difference is most apparent for the 
bleed boundary condition where the overall massflow 
through the region is specified in the PARC code while a 
velocity profile can be specified with the PROTEUS code. 

The standard turbulence model in the PARC code, just like 
the PROTEUS code, is the Baldwin-Lomax model. An 
alternative version of the code includes the two-equation K- 


e model which handles boundary layer history effects. For 
certain applications with flow separations, the more compli- 
cated K-e model has proven superior to algebraic models, 
reference 10. The K-e version of the PARC code is more 
computationally intensive per time-step. But with the better 
turbulence model, the solution may converge faster and 
thereby partially offset the slower time-steps. 

Another major difference between the codes deals with their 
time-accurate features. PROTEUS was conceived with 
time-accuracy in mind; while only recently were input 
options added to the PARC code to calculate time -accurate 
flows. This PARC option actually uses a ‘pseudo’ Runge- 
Kutta algorithm and is significantly slower than the standard 
Thomas algorithm. In either case, detailed comparison of 
unsteady flow solutions to experiment have not been 
reported. Work is actively being pursued to simulate an 
unsteady transonic diffuser flow with the PROTEUS code, 
however. In the present application, no attempt was made to 
use the time- accurate options of either code. Unsteady 
flows were simulated using the standard first-order time- 
accuracy of the codes with a globally constant time-step. 


III. Application - Supersonic Inlet Analysis 

This section describes both the initial preparation and the 
iterative process necessary to evolve a solution for the VDC 
inlet flowfield. First a discussion of the general procedure is 
given. More detailed discussions of the flowfield boundary 
conditions and the computational grid development are then 
presented. The care taken in selecting the boundary condi- 
tions and grid and the procedure used to evolve the solution 
has a dramatic influence on the quality and speed of the 
results. 


A. Procedure 

The study initially began by using the PROTEUS code with 
laminar viscosity. The initial flow field was set to Mach 2.5 
ffeestream conditions throughout the flow field and zero 
velocities at the inlet’s centerbody and cowl. The compres- 
sor face boundary is initially set as an extrapolation condi- 
tion. This setup should allow the inlet shocks to develop, the 
flow to compress nearly to critical conditions in the throat 
and then reaccelerate to supersonic conditions down 
through the diffuser and out the compressor face boundary. 
Once this flow solution reaches steady state conditions, var- 
ious levels of outflow “back” pressure are applied to posi- 
tion the normal shock downstream of the throat. 

This back-pressuring process is not straight forward. Since 
the change in back-pressure, (or any boundary condition 
change), occurs across some element of computational time, 
the change is an inherently unsteady event. Essentially, a 
change in pressure corresponds to an increase in momentum 
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due to the suddenness or acceleration of pressure change. If 
the pressure change occurs over a single iteration step as it 
does with the PARC code, a large transient shock forms 
whose strength is inversely proportional to the computa- 
tional time-step. This shock is analogous to an inlet ham- 
mershock that occurs in real supersonic inlet-engine 
systems when the engine stalls, references 20 to 22. Ham- 
mershock over-pressure, or transient pressures greater than 
freestream total pressure, develop behind the shock. If the 
transient is powerful enough, the inlet will unstart. Increas- 
ing the outflow boundary pressure must be done carefully 
and incrementally to avoid hammershock unstarts. 

A final note on procedure is needed to help define conver- 
gence to steady-state. Traditional levels of measured steady 
state flow convergence, (i.e. L2 residual reduction), were 
not reached in this application. This fact is due to the natural 
unsteadiness associated with strong shock turbulent bound- 
ary layer interactions. A global measure of convergence, 
mass flow conservation, was used in this study as well as 
other inlet studies, reference 16. Most of the computational 
time is required in positioning the terminal shock in the sub- 
sonic diffuser with the proper amount of back-pressure. The 
solution was examined periodically to check shock position 
and transient shock-induced oscillations. This examination 
often resulted in interactive adjustments of time-stepping, 
artificial viscosity and exit pressure. 

B. General Boundary Conditions 

For both the PROTEUS and PARC code the input boundary 
conditions were identical with the exception of the bleed 
slot boundary. The computational domain and boundary 
conditions are shown in figure 4. The flow was simulated 



FIGURE 4., Computational grid 

fine mesh, 447x199 points 

slightly ahead of the inlet centerbody spike, (the inflow 
plane), downstream to the location of an engine compressor 
face, (the outflow plane). The inflow plane is set to uniform 
freestream conditions. The outflow plane has an extrapola- 


tion boundary condition set initially. This allows the com- 
puted flow to reaccelerate to supersonic conditions through 
the outflow plane. The top and bottom contours are dictated 
by the centerbody and cowl internal wall surfaces. No slip 
wall boundary conditions are applied to these surfaces for 
the full Navier-Stokes solutions, while slip walls are applied 
for Euler solutions. The freestream boundary ahead of the 
cowl was simply set to the cowl lip radius and carried for- 
ward to the inflow plane. An extrapolation boundary condi- 
tion was applied here allowing the centerbody spike shocks 
to exit the computational domain with minimal reflection. 
For the axis of symmetry boundary that runs between the 
inflow plane and the centerbody spike, some difficulty with 
the PROTEUS code was encountered when the radius was 
set equal to zero. Consequently, a small diameter cylinder of 
1% of the cowl lip radius was artificially added to the geom- 
etry. The cylinder eliminates the R=0. line from the compu- 
tational domain and allowed the calculation to proceed. This 
geometry alteration was kept for the grid used with the 
PARC code. 

The inlet’s slotted bleed region was on the centerbody from 
x/R c =2.5831 to x/R c =2.6883. From a previous inlet bleed 
study, reference 15, the proper boundary condition for the 
bleed region was crucial to even obtain convergence in the 
solution. Of course, the accuracy in the flow solution was 
very sensitive to this boundary condition since it cut through 
flow with high pressure and velocity gradients. 

The implementation of bleed differs between the two codes 
used. In the PROTEUS code, individual grid point boundary 
conditions may be specified. The laminar bleed study done 
by Hamed, reference 19, gave profiles across the slot of tan- 
gential and normal velocity. This profile was used to specify 
the individual velocity components in the bleed slot region. 

For the PARC code a more restrictive boundary condition is 
applied over the entire bleed slot by setting a mass flow. The 
PARC code then sets a constant exit pressure over the entire 
slot length to deliver this mass flow. Errors in the calculated 
mass flow versus the input mass flow are used to adjust the 
exit pressure for the next iteration. This boundary condition 
is more restrictive than the one use in the PROTEUS code 
for the following reason. As a disturbance such as a shock 
traverses the bleed slot and causes a change in bleed mass 
flow, the PARC condition will try to maintain a constant 
mass flow while the PROTEUS condition will allow 
changes. This becomes important when the terminal shock 
is positioned by setting a compressor face pressure. 

C. Inlet Geometry/ Gridding 

To adequately resolve the various shock, boundary layer, 
separation, and bleed phenomena and their mutual interac- 
tions, a fairly fine mesh was required to analyze this inlet. 
Initially, the PROTEUS code’s grid size limitation was 
99x99. This was increased to 149x99 as a coarse mesh and 
to 447x199 as a fine mesh for this study. For the Proteus 
code, the coarse grid required about 3 million words, 
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(3Mwords), while the fine grid required nearly 20Mwords 
of memory storage. The PARC code has more efficiency 
array utilization and required only about 3Mwords of mem- 
ory for the fine grid. The coarse mesh cases were run on the 
LeRC Cray X-MP and Y-MP computers while the fine mesh 
cases were run on the Y-MP. With the coarse mesh, only 3 
points would be placed in the bleed slot region of the inlet 
throat with a uniform axial mesh. More grid was needed in 
this critical region to correctly model the complex flow phe- 
nomena. To accomplish this task, a simple two-dimensional 
algebraic grid generating program was written and used to 
adapt the mesh by packing and skewing control features. 

To make efficient use of the available number of grid points, 
shocks were resolved over a minimum physical space by 
skewing the grid along their expected positions. The 
remaining grid points were packed in the inlet throat and, in 
particular, the bleed region. 

To adapt the grid, the computational domain was sectioned 
into various regions. For example, figure 5 shows five 
regions: the freestream, the centerbody shocks, the pre- 
cowl, the cowl shock-throat-bleed, and the subsonic diffuser 
sections. Each of these regions have differing grid require- 



FIGURE 5., Blocking scheme to generate 
packedand skewed mesh. 

ments. The freestream and pre-cowl sections need minimal 
grid since the flow is uniform or nearly uniform within 
them. The centerbody shocks region needs moderate grid- 
ding packed near the inflow and outflow boundaries to 
resolve the shocks. The cowl shock-throat-bleed region 
needs maximum packing to resolve the shock/boundary 
layer/bleed interaction. Finally, the subsonic diffuser needs 
moderate gridding packed on the upstream boundary. 

The sectioning of the computational domain occurs by spec- 
ifying quadrilaterals. Each section must have a side coinci- 
dent with: an inflow and outflow crossplane, the cowl 
surface, and the centerbody surface. Grid lines can be set up 
along expected shock wave positions by skewing the inflow 
and outflow crossplanes. Control coefficients for the inflow 
cowl point and centerbody point determine packing within a 
section. The outflow packing is faired with the downstream 
sections inflow control coefficients. By packing the inflow 
and outflow crossplanes, the center of the section becomes 
sparse and vice versa. 

This scheme results in a distribution of crossplane end 
points defining the geometry of the computational domain, 
the packing and the skewness. Finally, to maintain an ade- 
quate degree of smoothness, (change in axial station <10% 


to 20% is a good rule of thumb), a simple central difference 
averaging scheme adjusts the ‘non-smooth’ control points 
iteratively until they all meet the smoothness criteria. Note 
that this procedure does not specify the radial distribution of 
points near the cowl and centerbody surfaces. Radially 
packing of grid points is done independently within the 
PROTEUS code. 

Figure 6 presents some results of the grid generation pro- 
cess. The initially packed and skewed mesh is shown in the 
above grid. This grid was smoothed so that the maximum 
change in streamwise packing of the points did not exceed 
10%. For the 10% smoothed mesh, calculations using the 
PROTEUS code intermittently failed. A final grid packing 
smoothness of 5% proved robust and was used in all subse- 
quent PROTEUS and PARC calculations. This 5% 
smoothed grid is shown in figure 4. The mesh shown is for 
the fine mesh. The coarse mesh was similarly developed 
using 5% smoothing. 



FIGURE 6., Smoothed mesh 

Using this grid generation tool, the number of points in the 
throat bleed slot increased from 3 to 7 for the coarse mesh, 
(149 axial points). For the fine mesh, (447 axial points), a 
total of 53 points were placed in the bleed slot. 


IV. Results 

Results of computational analysis for the VDC inlet are bro- 
ken into eight areas. The analysis will be compared to static 
pressure distributions and pitot rake measurements from the 
fixed hardware tests. In addition to the comparison, data 
from the computational analysis will be presented as Mach 
number and pressure contour plots. 

The results from the FNS analysis code will be discussed at 
length. The discussion of the FNS code predictions are 
divided into four areas: A) Euler and laminar cases, B) Tur- 
bulent cases with bleed, C) Design back-pressure case, 
D) Effect of cowl bleed. Along with these topics, the need 
and impact of cowl translation will be discussed. 

A. Euler Cases 

Results from the FNS analysis are based on the two similar 
codes, PROTEUS and PARC. Initial results were obtained 
just using the Euler subset of the FNS equations and using a 
fairly conventional mesh. The mesh for this case was uni- 
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formly distributed in the radial direction and slightly packed 
in the streamwise direction. The grid dimensions were 
99x99. This case had extrapolation outflow boundary condi- 
tions allowing the flow to reaccelerate to supersonic condi- 
tions downstream of the throat. 

Mach number contours are shown in figure 7 from this case. 
The results were obtained with the PROTEUS code. Exami- 
nation of the contours shows significant shock wave smear- 
ing in the physical domain. In fact, the cowl shock is not 
sharply defined, and the entire cowl compression appears to 
be distributed both well upstream and downstream of the 
shoulder. (Recall that the inlet design was for shock cancel- 
lation and focussed cowl compression at the centerbody 
shoulder). 



An effort in skewing and packing the grid as well as adding 
more grid points to the calculation domain was strongly 
motivated by the previous result. Figure 8 shows Mach 
number contours resulting from intermediate calculation of 
the inviscid flowfield with a refined mesh. Also shown is the 
previous calculation done on the uniform mesh. Although 
both of these calculations eventually unstarted, the sharp- 
ness of the shocks is clearly much better for the skewed, 
packed mesh. 



FIGURE 8., Effect of «grid placement on 
developing fculer solution 


Calculations continued for the coarse 5 % mesh for the 
design inlet geometry without bleed. Euler and laminar cal- 


culations failed since they predicted the cowl shock inter- 
section forward of the shoulder. This forward intersection 
caused a separation at the centerbody surface that grew and 
eventually caused the inlet flow to unstart. From the fixed 
hardware tests, at least a 2% bleed mass flow ratio was 
needed to maintain a started inlet. This result might help 
explain why the viscous calculated flowfield unstarted. The 
Euler result is more puzzling, since the prior MOC and PNS 
calculations suggest a started, shock-on-shoulder inlet flow 
should have been predicted. These calculations were done 
by the PROTEUS code and the coarse mesh size was lim- 
ited by the memory size of the available Cray-XMP super- 
computer. (The YMP machine with its larger memory was 
not available at that time). 

Restricted in mesh size, additional cases were run by first 
overspeeding the inlet and then by translating the center- 
body. By overspeeding the inlet, the higher Mach numbers 
both caused lower shock angles as well as a higher throat 
Mach number. With the lower shock angle the intent was to 
cause the shock to hit downstream of the shoulder. Unfortu- 
nately, the conical centerbody shocks were also more 
oblique and were ingested at the cowl lip plane. Figure 9 
shows Mach number contours for this case. The centerbody 

Mach 3.0, 'Oversped - Design Geometry 



Closeup of Throat Region 



FIGURE 9., M3 oversped laminar 

solution, Mach contours. 

shock ingestion caused separation on the cowl surface. The 
multiple shocks from this cowl separation do however inter- 
sect the centerbody at, and downstream of the shoulder. For 
this case the cowl separation grows and eventually cause 
inlet flow unstart. 

A similar M3 oversped case was run, but now with the cen- 
terbody translated forward 0.063R C relative to the cowl sur- 
face. The centerbody shocks were no longer ingested and 
the cowl boundary layer remained attached, (figure 10). 
Note also, the cowl shock hits well downstream of the 
shoulder. It passes through the strong expansion of the 
shoulder and causes the centerbody boundary layer to sepa- 
rate. This case did remain started 
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Mach 3.0. 'Oversped' - Cowl Translated 


B. 


Turbulent Case with Bleed 




FIGURE 10., M3 laminar solution with 
cowl translation. 

The next case attempted was at the design Mach number of 
2.5, cowl translated at 0.063R C and about 6% centerbody 
bleed mass flow ratio. At this time, the CRAY-YMP became 
available, and solution with the fine mesh was pursued. 
Since the exit massflow reached a constant value, this case 
reached a steady started inlet flow prediction. Mach number 
contours of figure 1 1 show the cowl shock hitting within the 
bleed slot, as well as a flow expansion caused by the bleed 
flow leaving the computational domain. Due to this expan- 
sion, the flow turns towards the centerbody and a fairly 
strong oblique shock forms as the flow hits the centerbody 
at the downstream edge of the bleed slot. These results were 
for the PROTEUS model and used the velocity profile from 
the Hamed and Lehnig study, reference 19, to simulate the 
bleed boundary. 


MACH 2.5. Cowl Translated Aft, 
Slot bleed -6% 



COMTOUN I tin I 



As mentioned, the started laminar viscosity case with 6% 
centerbody bleed and cowl translation did converge in terms 
of exit massflow. This case was initially chosen to explore 
the process of back-pressuring, or changing the outflow 
boundary condition. For the design geometry, figure 12 
shows Mach number contours in the subsonic diffuser for 
various computational time slices as the outflow boundary 
pressure is increased. The terminal shock advances forward 
towards the throat of the inlet. The exit static pressure ratio 
was set to 1 1 .62 and global time-stepping was used when 
the shock approached the throat. 



FIGURE 11, M2. 5 laminar solution 
with cowl translation. 


With these encouraging results, emphasis was shifted back 
to the design geometry without cowl translation. The grid 
was suitably modified and the earlier solution was used as 
initial conditions for the design geometry. Solutions were 
obtained for the turbulent flow case of the design geometry. 


FIGURE 12., Effect of back-pressure, 
subsonic Mach contours. 


Global time-stepping locks a constant time-step to advance 
the solution throughout the computational domain. Local 
time-stepping advances the solution for each grid point at 
different rates depending on the local value of CFL number 
which sets numerical stability. Local time-steps generally 
converged to a steady state solution faster than global time- 
stepping. However, this faster convergence seems to create 
larger oscillations in the solution when the outflow bound- 
ary condition is changed. When the terminal shock 
approached the throat, (critical inlet flow), even small pres- 
sure oscillations of 5% could push the shock upstream of 
the throat and cause the solution to unstart. The global time- 
stepping helped damp out spurious transients in the solution 
induced during local time-stepping. 

In setting the value of exit pressure level, guidance was 
gained from the experimental data. The data indicates the 
correct level of pressure ratio for the outflow condition 
ranges from: 13.7 for supercritical, 14.1 for critical, to 14.4 
for “peak” or subcritical inlet flow. This data was for about 
2% centerbody bleed flow and a slightly different subsonic 
diffuser. The diffuser had vortex generators to prevent sepa- 
ration installed downstream of the throat. Additional data 
was found in reference 2 that showed the presence of the 
vortex generators affected the outflow pressure ratio by only 
3%. This brings the above-mentioned range to: 13.3 for 
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supercritical, 13.7 for critical, and 14.0 for “peak” or sub- 
critical inlet flow. However, all of these levels of back-pres- 
sure led to unstarted flow computationally when they were 
suddenly imposed at the exit plane. A more conservative 
pressure ratio of 1 1 .62 was arrived at by trial and error. . 



FIGURE 13., Flow oscillation, Mach 
contours, Pcf/Pinf=11.62 


At this pressure ratio of 11.62, the solution did not reach 
any form of steady state solution. In examining the solution, 
the terminal shock’s final position was found to oscillate 
around a fixed location, figure 13. This solution was 
obtained using the PROTEUS code with turbulent viscosity. 
From the freestream entrance plane to nearly one throat 
height downstream of the aft edge of the centerbody bleed 
slot, the Mach number contours remain constant with 
respect to computational time. Just downstream, the termi- 
nal shock location first advances forward and then collapses 
back downstream. Figure 14 shows the centerbody static 
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solution 
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(447x199) 
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oscillation is 
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( 14 hrs Cray-YMP) 


FIGURE 14., Flow oscillation, Pcf/Pinf=11.62 
Centerbody Static pressure distribution. 


pressure distribution for these cases. Again, the shock loca- 
tion indicated by the sharp pressure rise around X/R c =3 first 
travels upstream, then downstream. Also, note the rise is 
much sharper for the upstream traveling shock compared to 


the downstream traveling shock. The period of this oscilla- 
tion took over 14,000 iterations or 14 CRAY-YMP Cpu 
hours to compute. Experimentally, dynamic distortion, a 
measure of flow unsteadiness, was found to increase at 
supercritical condition. This result supports the flow oscilla- 
tion phenomena encountered with the FNS analysis. 

Another point to mention on the centerbody pressure distri- 
bution, figure 14, is the pressure rise at X/R c ~2.4. This axial 
location is upstream of the shoulder, X/R c =2.5831, and indi- 
cates the cowl shock is impinging forward of the shoulder. 
The data also indicates this phenomena occurs. A more 
detailed view of this region of the flow is discussed next. 

The PARC code was run for the same case with the same 
grid. Comparison of Mach number contours is shown in fig- 
ure 15. The contours are for a detailed region through the 
inlet throat. Again, upstream of the terminal shock location, 
agreement between the two codes is good. Evident is the 
cowl shock hitting slightly forward of the shoulder causing 
a small separation. The cowl shock reflects from this separa- 
tion and then crosses back to the cowl surface. The reflec- 
tion on the cowl surface is of sufficient strength to separate 
the cowl boundary layer and cause a Mach reflection. 
Downstream of the Mach reflection’s normal shock, the 
flow is slightly subsonic but accelerates back to supersonic 
conditions. The PARC code resolves this phenomena more 
crisply than the PROTEUS code. The shock continues to 
reflect and coincides with the oblique shock at the aft end of 
the bleed slot on the centerbody surface. In the PROTEUS 
analysis, this oblique is weak and the flow remains super- 
sonic downstream. The PARC code predicts a strong 
oblique shock that coalesces with the terminal to generate 
subsonic flow. For both code predictions, the terminal shock 
is locally unsteady. Perhaps the most important difference 
between the code is the speed and memory requirements. 
For this case, PARC used only 40% of the Cpu time per iter- 
ation and 15% of the memory. 



To further display the features of this flow, figure 16 shows 
the particle traces for this case from the PARC code analy- 
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sis. The small separation forward of the centerbody shoul- 
der bleed slot and the bleed flow exiting through the slot are 
clearly evident . The result also suggests a need for addi- 
tional static pressure instrumentation to better quantify the 
extent of the separation. . This case was done with the cen- 



FIGURE 16., M2.5 design case, Particle 
traces, Pcf/Plnf=11.62 


terbody bleed flow rate at 2%, so its solution is directly 
relatable to the data. 

Figure 17 shows the results of the comparison of pitot pres- 
sures between the experiment and numerical solution. The 
various rakes are located on figure 14. Comparison is excel- 
lent for the cowl throat rake. The flow on the cowl to this 
point is unaffected by separations and by the centerbody 
bleed, so boundary layer growth and oblique shock pressure 
level should be correctly modeled by the code. Pitot pres- 
sure profiles for the other two rakes show moderate agree- 
ment. These rakes are downstream of separations and the 
bleed slot. Therefore, they are strongly affected by phenom- 
ena that are, at best, only approximately simulated by the 
turbulence model and the bleed boundary condition. 

VDC INLET, PARC vs. DATA 





FIGURE 17., Pitot pressure rake data 
versus prediction. 


C. Design Back-pressure Case 

The next few figures show the computations for an exit 
pressure ratio of 13.7 that predicts unstarted inlet flow. Fig- 
ure 18 shows the static pressure distributions as the outflow 
boundary pressure is suddenly increased. The terminal 
shock pressure-rise moves upstream with time. Downstream 
of the shock, static pressure levels oscillated in computa- 
tional time. 



FIGURE 18., Design back-pressure, 

static pressures, sudden change 

The plot of total pressure distribution is shown in figure 19. 
These values are for a line midway between the cowl and 
centerbody surfaces. The problems in the transient pressure 
pulse are evidenced by increases in total pressure ratio. 
(Note, total pressure is inversely proportional to entropy for 
a steady ideal gas. For a transient flow, however, entropy 
can decrease in the streamwise direction for a moving 
shock, see Appendix A). If a steady state solution is ulti- 
mately sought, any increase in total pressure can be con- 
strued as an indication that the flowfield is unsteady. Since 
both PROTEUS and PARC are solving the unsteady Navier- 
Stokes equation, albeit inaccurately in time, the pressure 
oscillation can occur and extreme care must be taken to pre- 
vent overwhelming pressure oscillations. 

The sequence in figure 20 shows the end result, the calcu- 
lated flowfield unstarts. A total variation diminishing 
scheme, (TVD), which reduces numerical oscillations about 
shock waves, (references 17 and 18), might be a more suit- 
able code to do inlet back-pressure studies. 

Because of the extreme sensitivity of shock position to the 
exit pressure boundary condition, further effort was made to 
gradually increase exit pressure and simulate critical inlet 
operation. PROTEUS was chosen for this effort due to two 
reasons. Since exit pressure can be changed gradually over a 
specified number of iterations, the problem of hammer- 
shock overpressure can be greatly reduced. The secomd rea- 
son for using PROTEUS lies with its more flexible bleed 
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FIGURE 21. .Design back-pressure, 

Mach contours, slow change. 


Mach number contours which do not exhibit a Mach reflec- 
tion are a consequence of the reduced internal contraction. 


The integrated massfiow was 97.5%± 0.2%. Total pressure 
recovery was 88% ± 0.2% and the bleed massfiow rate was 
-1.9%. Also shown in figure 21 is the comparison of static 
pressure distribution. The agreement with data is very good, 
the slight descrepencies are due to greater calculated loss in 
total pressure recovery, (the experimental recovery was 
89.5%). 


FIGURE 20., Design back-pressure, Unstart 
Mach contours, sudden change. 

model. Both of these reasons were crucial to the successful- 
simulation of critical inlet operation, figure 21. 

For this case the exit pressure was increased from 12.74 to 
13.7 over 20,000 iterations. Even with this slow change of 
pressure a slight overpressure occurs, travels upstream, and 
moves the terminal shock forward onto the bleed slot The 
bleed boundary condition model, which was based on 
Hamed’s velocity profile, did allow the massfiow to increase 
by -50% during this terminal shock/bleed interaction. The 
shock eventually moved back downstream to the aft edge of 
the bleed slot and remained stable. 

Also, this case was run with slight amount of cowl transla- 
tion, X/R c =0.035. This translation was picked to position 
the cowl shock intersection with the centerbody down- 
stream onto the bleed slot. The intent was to eliminate the 
small separation bubble upstream of the slot. The result, 
shown in figure 21, was a decrease in the extent of the sepa- 
ration. Because of the translation, approximately 3% more 
flow is spilled and the internal contraction is reduced. The 


D. Effect of Cowl Bleed 

In order to control the oscillation of the terminal shock, a 
region of cowl bleed was simulated with the PARC code by 
specifying mass flux through the cowl boundary. The VDC 
inlet’s porous cowl bleed is located between x/R c =2.556 to 
x/R c =3.105. An initial run was attempted with the entire 
area controlled by the bleed boundary condition. 

Unfortunately, this large, semi-free boundary seemed to 
aggravate rather than reduce the flow unsteadiness, (figure 
22). In fact, the momentum vectors at the cowl bleed bound- 
ary flucuate wildly both spatially and temporally. The vec- 
tors for the case presented indicate alternating areas of 
blowing and bleeding through the ‘bleed’ boundary. In the 
Mach number contours, a very large separation has devel- 
oped on the cowl surface due to the blowing momentum at 
the end of the region. 

This result suggests the need for a more refined simulation 
of the bleed region. Limiting the extent of the bleed region 
may help direct the flow and reduce the areas of blowing. 
Compartmentalizing the bleed, or having two independently 
specified bleed regions, is also worthy of investigation. Fur- 
ther work in refining the cowl bleed is recommended, par- 
ticularly because of its potentially dramatic impact on the 
upcoming wind tunnel testing. 
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FIGURE 22., Effect of cowl bleed, 

entire bleed region open. 


V. Conclusions 

Summary: 

A series of analytical tools were applied to understand the 
flowfield in a supersonic inlet. These tools had varying 
degrees of success when compared to existing data from a 
wind tunnel test of the inlet geometry. Some of the major 
results are summarized below: 

1) Elliptic Euler analysis failed due to instability 

in the centerbody flow separation, which led 
to an unstarted flowfield. 

2) FNS viscous analysis predicted unstart for the 

design geometry without centerbody bleed, 
agreeing with experimental testing. Both 
analysis and testing agree that the design ge- 
ometry will remain started at a bleed mass 
flow ratio, m bl /m c , of 2%. Analysis further 
predicts that the inlet will remain started 
without bleed at a centerbody translation of 
X/R c =0.065. 

3) For a supercritical inlet flowfield, FNS turbu- 

lent viscous analysis with bleed agreed well 
with experimental data forward of the cen- 
terbody bleed slot including a region of sep- 
arated flow just upstream of the slot. 

Downstream of these regions the agreement 
is fair, pointing out deficiencies in the bleed 
and turbulence modeling. These deficiencies 
prevented analysis of critical and peak inlet 
conditions. 

4) FNS codes indicated an inherent unsteady flow 

oscillation of the terminal shock at a super- 
critical flow condition. This agrees with em- 
pirical experience that suggest that dynamic 


distortion increases as the terminal shock 
moves downstream from the inlet throat. 

5) Both FNS codes failed to correctly predict the 

inlet recovery. Analytical predictions were 
-1.5% below measured recoveries. 

6) Preliminary analysis with the PARC code indi- 

cated that a fully open cowl bleed pattern 
may enhance rather than eliminate flow un- 
steadiness and subsonic diffuser separations. 

7) Analysis predicted that the flow separation up- 

stream of the centerbody shoulder can be 
virtually eliminated with a centerbody trans- 
lation of x/R c =0.035. 
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The analysis for a one-dimensional transient duct flow 
shows that peak pressure during a transient can be substan- 
tially greater than the steady-state total pressure of the flow. 
This result contrasts to the three-dimensional steady-state 
result: total pressure can never increase in the flow direction 
provided no work is added to or heat is taken from the flow. 

A schematic of the flow phenomena is given below, figure 
A.l 



FIGURE A.1 , 1 D Transient duct flow 


The assumptions made in this analysis are one-dimensional 
ideal air flow with y=1.4. Also assumed was that zero 
velocity occurs after the transient and that a normal shock 
travels upstream at a velocity, u s . The assumption of zero 
velocity ahead of the blockage valve, (or stalled engine), is a 
reasonable simplification, but it does not represent the worst 
case scenario. Flight experience has proven that during a 
hammershock transient, reverse flow can occur in the pro- 
pulsion system, i.e. flames shoot out of the inlet! For this 
severe situation, over-pressures beyond that indicated by the 
following simple analysis can occur. 

The initial frame of reference is fixed to the duct walls. 
Now, if we change our reference frame to traveling with the 
moving normal shock, the steady-state normal shock rela- 
tions can be applied. In the new frame of reference, the 
Mach number ahead of the shock becomes: 


The Mach number ahead of the shock in the wall fixed 
frame, M a , is related to the shock fixed frame Mach number 
by: 


JW a 


M 


a 



VII. Appendix 

Appendix A. Hammershock Over-pressure Analysis 

Equations are developed for the one-dimensional transient 
duct flow undergoing a sudden blockage at the exit. This 
analysis demonstrates the concept of hammershock over- 
pressure. 


The pressure ratio, downstream static to upstream total, is 
thus: 

Also, the Mach number downstream of the shock in the wall 
fixed frame, M s , is related to the shock fixed frame Mach 
number by: 
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Appendix B. Axisymmetric Governing Equations 

The Full Navier-Stokes equations for axisymmetric flow are 
given along with assumptions used in analyzing the VDC 
Inlet. Also presented are the equations for the Baldwin- 
Lomax algebraic turbulence model. 

The analysis used in the FNS codes for the VDC Inlet study 
are the well known Navier-Stokes equations written in 
cylindrical coordinates for axisymmetric flow. These equa- 
tions are somewhat more complicated than the classic two- 
dimensional planar equations due to the dependence of vari- 
ous terms on the radial direction. In vector notation the axi- 
symmetric equations with swirl become: 

f , «® +Tr (,f> + " * £i friv) 


If we use M as as a parametric variable, we can solve the 
equations for pressure ratio and M a and plot the resulting 
relationship, (figure A.2) 

Pressure Ratio, P s / P t a 

A 



FIGURE A.2, Hammershock 

Over-Pressure Ratio 

Note that the maximum over-pressure reaches a level of 
1.84 times the freestream total pressure at slightly less than 
sonic conditions in the duct, figure A.2. As the Mach num- 
ber further increases, the over-pressure begins to decrease as 
the normal shock losses overwhelm the transient increase in 
total pressure. 


Unsteady terms: 


Q = [p pu pv pw E ~ t ] 


Axial derivatives: 


P u 
pu 2 + p 
p uv 
P MW 

(E T + p)u 


Radial derivatives: 


F 


pv 
p uv 

pu 2 + p 
pviv 

(E T +p) v 


A typical supersonic inlet has nearly sonic conditions in the 

throat. This analysis clearly shows how over-pressures Swirl derivatives: 

greatly exceeding the freestream total may develop during 

transient inlet phenomena. Therefore, when an engine stalls 

downstream of an inlet, or when the compressor face 

boundary condition is rapidly changed in a computational 

study, hammershock over-pressures are produced. 
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0 
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-(P ™ 2 + p) 
pvw 
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Axial viscous derivatives: 


Ey ~ IT 


xQ 


UX + VT +WT <7, 

xx xr xQ p r 


Radial viscous derivatives: 


F * = k 


re 


UX + VT +WX Q -— q, 
xr rr rQ p r ^ r 


,3w W 


i, - 




In the equations, x, r, and 8 are the axial, radial and circum- 
ferential coordinate directions, respectively. The velocity 
components: u, v and w, are along those respective direc- 
tions. The equations have been nondimensionalized by ref- 
erence quantities of length, velocity, density, temperature, 
viscosity, and thermal conductivity. Pressure and total 
energy are 

L 

divided by p u , and time by — 

r r u 


The reference Reynolds and Prandtl numbers are thus 
defined as: 


Swirl viscous derivatives: 


Re, = 


p,“A 

Pr 


H v = 



0 

0 



0 


Due to the extra complexity introduced by the axisymmetric 
terms, the strong conservation law form of the equations is 
lost. The viscous and heat transfer terms are: 


x 


XX 


~ du . 
2ii— + X 
dx 


~du 1 3(rv) 
r dr 


>] 


Pr 


r 


Pr“ 


2 

r 


kjt 


The above equations, which are written in Cartesian coordi- 
nates, are then transformed to a general coordinate system. 
Finite difference techniques are then applied to discretize 
these equations. The Beam and Warming method is used to 
difference the equations in time. For both the PROTEUS 
and PARC codes, first order accurate time differencing was 
used. Also, the spatial derivatives were centrally differ- 
enced. Details of the implementation of the differencing 
schemes of the two codes are given in references 8, 9 & 10. 


. rdu l 3 (rv) -| 


.v , r8« 1 d(rv) -] 


,du dv. 


% 


The equations are valid for laminar or Reynolds time-aver- 
aged turbulent flows. For the viscous terms, (I, X and k 
represent effective coefficients for turbulent shear flow. The 
Baldwin-Lomax turbulence model is used by both codes 
and, for wall-bounded flows, is given by: 

^tot = »lam + ^urb 


turb 


( Pr ) ou u^n>y b 
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for free shearing turbulent flows, **. = . 

yn is the distance normal to the wall for bounded flows 
or the minimum distance to the local extrema of velocity 
magnitude for free shear flows, y ^ is the minimum value of 
y at which ^ inner equals outer . For the 
PROTEUS code an averaging procedure is used if the shear- 
ing flows associated with the two walls merge. The PARC 
code uses the standard implementation strategy. 


I = y * /A ' ) 

where K is the Von Karman constant usually taken as 0.4. 
Also to find the effective thermal coefficient, the Reynolds 
analogy and a turbulent Prandtl number, Pr lurb , of 0.9 was 
used: 

. _ C p^turb 


The standard Baldwin-Lomax outer region is given by: 

^tKuter = KC cpP F Kleb F wake 

where K=0.0168, 0^=1. 6, and F wake is given by 


Finally, the second coefficient of viscosity is given by: 


X 


turb 


2 

3^(«r6 


F waU = min 


y F 

J max rr 


' max max 

2 

C wk v diffp 
^ max' 


— / 2 2 2 
V diff ~ (u + v 4- w ) - min 


J (u 2 + v 2 + w 2 ) 


F ( j y /A ) for wall-bounded flow \ 

V yjft| for free shear flow ) 

F max is the maximum of F(y n ) and y max is the value of y n at 
this position. 


In the analysis of the VDC Inlet the swirl velocity, w, was 
assumed to be zero and the total enthalpy was assumed to be 
constant. Therefore the circumferential and energy equa- 
tions drop out of the coupled matrix of nonlinear partial dif- 
ferential equations, considerably reducing the complexity of 
the computation. Also assumed for the equation of state is 
that the airflow behaves as an ideal gas with the ratio of spe- 
cific heats, y, equal to 1.4. Also Sutherland’s viscosity rela- 
tion for air is used: 


P"lam 


_ jlSl 


d+T suth ) 

(T+T suth ) 




/ dw vv 2 3w 2 3v du 
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where T sulh is the nondimensional Sutherland temperature. 
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and C wk =0.25, A + =26.0, B=5.5, 0^=03 


The standard Baldwin-Lomax inner region is given by: 
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